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Abstract 

In the CP-violating Minimal Supersymmetric Standard Model, we study 
the pair production of neutralinos at center-of-mass energies around the heavy 
neutral Higgs boson resonances. For longitudinally polarized muon beams, 
we analyze CP asymmetries which are sensitive to the interference of the two 
heavy neutral Higgs bosons. Due to radiatively induced scalar-pseudoscalar 
transitions, the CP asymmetries can be strongly enhanced when the reso- 
nances are nearly degenerate, as in the Higgs decoupling limit. Fhe Higgs 
couplings to the neutralino sector can then be analyzed in the presence of 
CP violating phases. We present a detailed numerical analysis of the cross 
sections, neutralino branching ratios, and the CP observables. We find that 
radiatively induced CP violation in the Higgs sector leads to sizable CP- 
asymmetries, which are accessible in future measurements at a muon collider. 
However, we expect that our proposed method should be applicable to other 
processes with nearly degenerate scalar resonances, even at hadron colliders. 
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1 Introduction 



The CP-conserving Minimal Supersymmetric Standard Model (MSSM) contains 
three neutral Higgs bosons [1-5], the lighter and heavier CP-even scalars h and 
H, respectively, and the CP-odd pseudoscalar A. While the MSSM Higgs sector is 
CP-conserving at Born level even in the presence of CP-phases, loop effects medi- 
ated dominantly by third generation squarks can generate significant CP-violating 
scalar-pseudoscalar transitions, leading to mixing of the neutral Higgs states into 
the mass eigenstates Hi, H2, H3, with no definite CP parities [6,8-14]. 

It is well known that mixing of states with equal conserved quantum numbers 
is strongly enhanced when these states are nearly degenerate, i.e., when their mass 
difference is of the order of their widths [9,10]. This degeneracy occurs naturally 
in the Higgs decoupling limit of the MSSM, where the lightest Higgs boson has 
Standard Model-like couplings and decouples from the significantly heavier Higgs 
bosons [15]. In the decoupling limit, a resonance enhanced mixing of the states H 
and A can occur, which may result in nearly maximal CP- violating effects [8-10]. 
The general formalism for mass mixing in extended Higgs sectors with explicit CP 
violation is well developed [11-14], and sophisticated computer codes are available 
for numerical calculations [16-18]. Detailed investigations of the fundamental prop- 
erties of the Higgs bosons, both phenomenological and experimental, will be crucial 
for the understanding of the mechanism of electroweak symmetry breaking. 

In previous studies of the CP-conserving [19-27] and CP-violating Higgs sec- 
tor [28-35], it was shown that the CP-properties and couplings of the heavy neu- 
tral Higgs bosons can be ideally tested in collisions. Such a muon collider 
is a superb machine for measuring the neutral Higgs masses, widths, and cou- 
plings with high precision, since the Higgs bosons are resonantly produced in the 
s-channel [36-38]. The well controllable beam energy allows the study of the center- 
of-mass energy dependence of observables around the Higgs resonances. In partic- 
ular, the beam polarization plays an essential role for analyzing the CP nature of 
the Higgs sector itself. Not only backgrounds can be reduced, but the CP-even and 
CP-odd contributions of the interfering Higgs resonances to the observables can be 
ideally studied, if the beam polarizations are properly adjusted [23-33]. 

Besides the initial beam polarization, the final fermion polarizations are essen- 
tial to probe the Higgs interference. The secondary decays of the final fermions 
enable their spin analysis, and additional final state polarization observables allow 
for a complete determination of the CP-properties of the Higgs bosons [25-28]. For 
final state SM fermions //, with / = r, b, t, such polarization observables have 
been classified according to their CP and CPTlj transformation properties [28]. For 
the production of neutralinos [26] and charginos [27] with longitudinally polarized 
beams, it has been shown that asymmetries in the energy distributions of their de- 
cay products are sensitive to the Higgs interference of the CP-even and CP-odd 

T is the naive time reversal t — > —t, which inverts momenta and spins without exchanging 
initial and final particles. 
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Figure 1: Resonant Higgs exchange in neutralino pair production. 



production amplitudes. Thus the couplings of the Higgs boson to the neutralino 
and chargino sector can be analyzed in the CP-conserving MSSM. 

In this work, we extend the study of neutralino production at the muon col- 
lider [26] to the CP-violating case. For neutralino production — > XiX^i with 
longitudinally polarized muon beams, we define a CP-odd polarization asymmetry. 
We analyze the longitudinal polarizations of the produced neutralinos by their sub- 
sequent leptonic two-body decays x°j ~~ * ££n, £ = e, fi, r, with n = R, L for £ = e, /i, 
and n = 1, 2 for £ = r. With the energy distribution of the leptons we can define 
a CP-even and a CP-odd polarization asymmetry for the neutralino decay, which 
probe the neutralino polarization. First results for neutralino and also chargino 
production in the MSSM with explicit CP violation in the Higgs sector have been 
reported in Ref. [7]. 

In Section [21 we give our formalism for neutralino production and decay with 
longitudinally polarized beams. In an effective Born-improved approach, we include 
the leading self energy corrections into the Higgs couplings. We give analytical 
formulas for the production and decay cross sections and distributions, and show 
that the energy distribution of the neutralino decay products depends sensitively on 
the Higgs interference. In Section [31 we define CP-odd and CP-even asymmetries of 
the production cross section and of the energy distributions. These observables are 
sensitive to the CP-phases in the Higgs sector, as well as to absorptive contributions 
from the Higgs boson propagators. In Section [H we present a detailed numerical 
analysis of the cross sections, neutralino branching ratios and the CP-observables. 
We give special attention to the \fs dependence of the observables and analyze their 
dependence on the CP violating phase 0^, of the common trilinear scalar coupling 
parameter A t = Ab = A T = A, and on the gaugino and higgsino mass parameters 
M.2 and /i, which mainly determine the Higgs couplings to the neutralinos. We 
summarize and conclude in Section [5) 



2 Neutralino production and decay formalism 

We study CP violation in the Higgs sector in pair production of neutralinos 

+ (1) 

with longitudinally polarized muon beams. The Feynman diagram for Higgs boson 
exchange is shown in Fig. [U We will analyze the process at center-of-mass energies 
of the nearly mass degenerate heavy neutral Higgs bosons H 2 and H 3 . They will be 
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Figure 2: Feynman diagrams for non-resonant neutralino production n + jj, — > 



resonantly produced in the s-channel. The significantly lighter Higgs boson H\ and 
the Z-boson are also exchanged in the s-channel, however far from their resonances. 
Together with the smuon exchange £il,r m the t- and w-channels, they contribute the 
non-resonant continuum to the neutralino production, see their Feynman diagrams 
in Fig. El 

To analyze the longitudinal polarizations of the produced neutralinos, we con- 
sider their subsequent CP-conserving but P-violating leptonic two-body decays 

+ e = e,»,T, (2) 

with n = R, L for £ = e or t = //, and n = 1, 2 for i = r. 



2.1 Lagrangians and amplitudes for Higgs exchange 

CP violation of the MSSM Higgs sector is induced by scalar-pseudoscalar mixing 
at the loop level. We will include these mixings effectively in the interaction La- 
grangians for neutralino production (pQ) via Higgs exchange Hj~, with k = 1,2, 3, 

= Ji[c H L k ^P L + c H R ^P R }iiH k , (3) 

^CCB = lu4™*PL + 4™*PdXjH h . (4) 

We obtain the effective Higgs couplings to the initial muons, c L k ^, and the final 
neutralinos, c^ k ^ X3 , from their tree level couplings 

C L,R ~ ^ka^L,R ' \°J 

c L*R Xj = C"ka c L,R Xj i h a = h,H,A, (6) 

with 

C = C~ lT . (7) 

The tree level couplings c^"^ and c^"^ Xj are defined and discussed, e.g., in Refs. [3, 
26]. The matrix C diagonalizes the Higgs mass matrix M at fixed p 2 . In the tree- 
level basis of the CP-eigenstates h, H, A, the symmetric and complex mass matrix 



4 



at momentum squared p 2 is given by [16] 



M(p 2 



( 



\ 



mi 



-Z hh (p 2 
-t hH {p 2 ) 
-t hA {p 2 ) 



m 2 H - t HH (p 2 ) 



-t hA {p 2 ) 
-t HA {p 2 ) 



-t HA {p 2 ) 



t AA (p 2 ) ) 



Here £ rs (p 2 ) with r,s = h,H,A are the renormalized self energies of the Higgs 
bosons at one loop, supplemented with higher-order contributions, see Ref. [16]. 
When the Higgs bosons are nearly mass degenerate, these corrections are enhanced 
by the Higgs mixing. The propagator matrix 



Ks( P 2 ) = -i[p 2 - M(p 2 )} 



1 

rs ' 



(9) 



has complex poles at p 2 = M. 2 Hk = M# —iM Hk T Hk , k = 1,2, 3, where M Hk and Y Hk 
are the mass and width of the Higgs boson mass eigenstate H k , respectively. We 
evaluate the mass matrix M(p 2 ) at fixed p 2 = Mfj 2 in its Weisskopf-Wigner form 
with the program FeynHiggs 2.5.1 [16, 17], in order to obtain the diagonalization 
matrix C, as well as the Higgs masses and widths. Here we neglect the momentum 
dependence of M(p 2 ), since, unless Higgs decay thresholds open at this energy, this 
dependence is weak in the resonance region p 2 w M]j 2 , . This approach corre- 
sponds to an improved-Born approximation. It includes the leading-order radiative 
corrections into the matrix C, but not the specific vertex and box corrections, as 
well as the subleading muon and neutralino self energy corrections. We give further 
details in Appendix [XI 

With the Born-improved effective couplings we write the amplitudes for neutra- 
lino production via Higgs exchange, see Fig. [Tj 



T p = A(H k ) [«(p M+ ) (c H L ^P L + c H R ^P R ) u{p,- 



x 



HP X o){cf kX ^P L + C H R ^P R )v(p x o] 



JlkXiXj 



(10) 



with the Breit-Wigner propagator for the Higgs boson 



A{H k ) 



11) 



s-M 2 Hk +iM Hk T Hk 

The Lagrangians and amplitudes for Z and £il,r exchange are given in Appendix [El 
and the Lagrangians for the leptonic neutralino decays (j2J) are given in Appendix O 



2.2 Squared amplitude 

In order to calculate the squared amplitude for neutralino production (CQ) and de- 
cay ([2]), we use the spin density matrix formalism of [40,41]. Following the detailed 
steps in Appendix [Dl where we also give the production amplitudes, the squared 
amplitude in this formalism can be written as 

3 

|T| 2 = 2|A(x?)| 2 (P- J D + $>?,£«), (12) 

a=l 
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with the propagator A(x°) of the decaying neutralino, see Eq. (ID. 291) . Here P de- 
notes the unpolarized production of the neutralinos and D the unpolarized decay. 
Ep and E£> are the corresponding polarized terms, and their product in Eq. (112j) 
describes the neutralino spin correlations between production and decay. With our 
definition of the spin density production matrix, Eq. (1D.30I) . Sp/P is the longitudi- 
nal polarization of Xp ^p/P is t ne transverse polarization in the production plane, 
and Ep/P is the polarization perpendicular to the production plane. We give ex- 
plicit expressions for the production terms P and E p in the next Section. The terms 
D and EJ, for neutralino decay are given in Appendix IFl 



2.3 Resonant contributions from Higgs exchange 

The expansion coefficients of the squared neutralino production amplitude ( Tl2|) sub- 
divide into contributions from the Higgs resonances (res) and the continuum (cont), 
respectively, 

P = P-cs + Pcont) Ep = E rcs + E cont , a = 1, 2, 3. (13) 

The continuum contributions P CO nt, E" ont are those from the non- resonant Z and 
P>l,r exchange channels, with P CO nt given in Appendix [E] and S" ont in [41]. 

In order to analyze the dependence of the resonant contributions P res and E^ cs 
on the longitudinal /i + and \i~ beam polarizations V + and P_ , respectively, we can 
expand^] 

P res = {l + P + V-)a + {P+ + P-)a 1 , (14) 
SL = (l + V + V-)b + (V+ + V-)b 1 . (15) 

Such an expansion proofs to be useful for discussing CP properties. The coefficients 
a and bi are CP-even, whereas a\ and b Q are CP-odd and vanish in the case of CP 
conservation [26]. They are given by 



a n = ^(2-fe)«n, b n = £>-Wn5 n = 0,l; 

k<l k<l 

with the sum over the contributions from the Higgs bosons Hk, Hi with k, I 
respectively, and 

s 



(16) 



1,2,3, 



a 



ki 



A 



(ki)\ 



\ c +\\c+\f ijC os(6; + S+ + 6 A ) 



— \c,, 



^RL 



mifrij cos(<5+ + 5. 



RL 

x 



(ki) 



ki 



-IA 



s 

V 

s 



(ki)\ 



n.— 1 1 pR-L 
1 1 X 



in 



iirij cos(5 + 5^ L + 5 A 



(*Q 



A(fc/) 


1 C~^~ \\c \ -\ 




cos(5+ - 




A(fc/) 




/Ay 


cos (5; - 





5a) 

5a ^ 



J (W) 

. (ki)' 



(17) 

(18) 
(19) 
(20) 



2 The resonant contributions Sj es and S 2 es to the transverse polarizations of the neutralino 
vanish for scalar Higgs bosons exchange in the s-channcl. 
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We have defined the products of couplings, see Eqs. 
neutralino indices % and j, 



and (jSj), suppressing the 



-'a(fcZ) 



H^ota Hiaa* 
Z R C R 

HkXXHixx* 



R 



C 



+ C 



HkXX^iXX* 



c£|exp(i<£; 



c 



R 



:f|exp(^f 



(ki) 
{ki) 



a 



the product of the Higgs boson propagators ffTTT) . 



KM) 



A(i7 fe )A(^; 



(ki) 



(21) 
(22) 

(23) 



and the kinematical functions fij = (s — m 2 — m 2 )/2 and Ajj, see Eq. (ID. 351) . 

Note that the coefficients bo and &i, which parametrize the neutralino polarization 
dependence in Eq. (1151) . vanish at threshold as \f\j- The coefficients a\ and b\ 
contribute only for longitudinally polarized muon beams. These coefficients are 
products of the Higgs boson couplings to the muons and neutralinos. Our aim is 
to determine these coefficients using polarization asymmetries and the production 
cross section. Since a muon collider provides a good beam energy resolution it is 
the ideal tool to analyze their strong y/s dependence. 

We neglect interferences of the chirality violating Higgs exchange amplitudes 
with the chirality conserving continuum amplitudes, which are of order m^/s/s. 
Note that contributions from Hi exchange will be small far from its resonance. 



2.4 Cross sections 

We obtain cross sections and distributions by integrating the amplitude squared 
|T| 2 (fT2"l) over the Lorentz invariant phase space element dPS 

da = —\T\ 2 dPS. (24) 

2(9 

We use the narrow width approximation for the propagator of the decaying neu- 
tralino. Explicit formulas of the phase space for neutralino production (CQ) and 
decay ([2]), can be found, e.g., in [42]. The /i + /i~-spin averaged cross section for XiX*j 
neutralino production is 



with the triangle function A^- f1D.35[) . and the average over the neutralino production 
angles in the center-of-mass system, 

P = — [ PdQ Y o. (26) 

47T J X i 

The integrated cross section for neutralino production ([T]) and subsequent leptonic 
decay x°j ~~ > ^ *n (0K with n = R, L for £ = e, /i, and n = 1, 2 for I — r, is given by 

* = (TW <^ {] ^^ P ■ D — &ij x BRtt? ,27) 

Explicit expressions for D are given in Appendix [F] Note that the integrated cross 
section <r" is independent of the neutralino polarizations. 
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2.5 Lepton energy distribution 



The differential cross section da (124j) . and thus the energy distribution of the lepton 
from the neutralino decay ([2]), depends on the longitudinal neutralino polarization. 
In the center-of-mass system, the kinematical limits of the energy of the decay lepton 
t = e, /i, r are 



with 



E max(min) = ^ e ± A e , (21 



2 2 

praax i ZTTmn 771 n — 777.? 

2 2m 2 o x j 



,2 ,^2 



pmax pmin 771 n — Tfl 



2m 2 o 



-|p x o|, £ = e,(x,T. (30) 



Using the definition of the cross section ( Tl2l) . ( f27|) . and the explicit form of T, 3 D ( IF. 541) . 
the energy distribution of the decay lepton ^ is [26, 27] 



do - 



e ± _ (J 



dE P 2A f 



Y? P [Ei — Ei 

1 + Tfci/tf 



P A, 



(31) 



with 7^± = =1=1. The factor 77" is a measure of parity violation in the neutralino 
decay. It is maximal r/^ = +1 and 77^ = —1 for the decay into en, fin and ^l^Li 
respectively, since the sleptons of the first two generations couple either purely left 
or right handed, if mixing is neglected. For the decay x? — > T± ^\2i the factor 



\hr |2 _ \ a f 12 



- \ b f 12, i f .12' w 

is generally smaller \rft\ < 1 due to stau mixing. The right and left Xj^n r couplings 
a T n j and are defined in (1G.60I) . 

Further in Eq. ( 13TI) . the coefficients P and £p of the squared neutralino produc- 
tion amplitude ffT2"]) are averaged over the neutralino production solid angle, denoted 
by a bar in our notation ([2"oT) . Due to the Majorana character of the neutralinos, 
the continuum contribution T? cont ( fT3l) is forward-backward antisymmetric [43], and 
vanishes if integrated over the neutralino solid angle. However, the resonant contri- 
bution E^ es (TlBl) from Higgs exchange is isotropic, thus 

^ 3p = An I ^r d ^x° = S res- ( 33 ) 



In Fig. [31 we show the energy distributions f[3"Tj) of the leptons from the 
decays — ¥ ^ + ^r an d X°j ~^ £~^r> for £ = e or /x The cutoffs in the energy 
distributions of the primary leptons £ + and i~ correspond to their kinematical limits, 
as given in Eq. ( 1281) . We see the linear dependence of the distributions on the lepton 
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Figure 3: Normalized energy distributions of the primary and secondary leptons 
for neutralino production [i + fi~ —> X1X2 an d subsequent decay chains X2 ~~ * ^~^~r-> ~^ 
Xiit (dashed) and x 2 -» ~ l R ~> X% ( solid ) for £ = e or /i. Here Af = -15%, 

corresponding to the masses of Table V+ = V- = 0.3, and ^/s = 494 GeV. 



energy. The slope of these distributions is proportional to the longitudinal neu- 
tralino polarization, and is solely due to the resonant Higgs contributions £j? cs . In 
addition, we show in Fig. [3] the energy distributions from the secondary leptons if 
from the subsequent decays — > I2X1 an d ^ ^Xi- 

Note that generally the parity conserving neutralino decays into the Z and the 
lightest Higgs boson, x° — > Z\\, and Xj ~^ HiXk' respectively, cannot be used to 
analyze the neutralino polarization. The resulting energy distributions are flat, due 
to the Majorana properties of the neutralinos. They imply that the left and right 
couplings obey \0'£\ = \0'£\ (EH. and \cf XjXk \ = \cf lXjXk \ ®. The decays 
are thus parity conserving, and therefore the parity violating factors analogous to 
77™ (IH^j) vanish. 



3 Asymmetries for neutralino production and de- 
cay 

In Eqs. (I14jl and (fl~5l) . we have expressed the resonant contributions P ms and S^ cs 
to the spin density matrix elements for neutralino pair production in terms of the 
longitudinal muon beam polarizations. In order to experimentally determine the 
four different combinations of products of couplings a , a\, b , and b\, we define 
asymmetries of the neutralino production cross section, as well as asymmetries of the 
energy distributions of the decay leptons. Together with the neutralino production 
cross section, these coefficients can then be experimentally determined, and thus the 
Higgs couplings to muons and neutralinos. 
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3.1 Asymmetry of the neutralino production cross section 



For the cross section of neutralino pair production = a{ji + ji — > XiX?)> Eq. ( |25l) . 
we define, for equal muon beam polarizations V + = V- = V, the CP-odd asymme- 
try [7,38] 

^■prod * s sensitive to the CP phases of the Higgs boson couplings to the neutralinos 
and to the muons. Denoting by T the naive time reversal t — > —t, which inverts 
momenta and spins without exchanging initial and final particles, the asymmetry 
^■prod i s a ^ so CPT-odd. Thus the asymmetry is due to the interference of the CP 
phases with the absorptive phases from the transition amplitudes. The absorptive 
phases are also called strong phases, and originate from intermediate particles which 
go on-shell. The asymmetry *4p° od is therefore sensitive to the CP phases of the Higgs 
boson couplings, as well as to the phases of the Higgs propagators. 

In the Higgs decoupling limit [15], the heavy neutral Higgs bosons are nearly 
mass degenerate. Thus a mixing of the CP-even and CP-odd Higgs states H and 
A can be resonantly enhanced, and large CP-violating Higgs couplings can be ob- 
tained [8,9]. In addition, CP phases in the Higgs sector lead to a larger splitting 
of the mass eigenstates H 2 and if 3. This in general tends to increase the phase 
difference between the Higgs propagators, giving rise to larger absorptive phases in 
the transition amplitudes. 



Using the definitions of the neutralino production cross section flzoj) . and of 
the P term 014p . we obtain 

jp°! - 2 J^1 (351 

Pr ° d " (1+P>0 + Pco„t' ^ ] 

We can thus employ the asymmetry to determine the CP-odd coefficient a\. The 
maximum absolute value of the asymmetry depends on the beam polarization V 

/lP°l,max _ ^' (op\ 
"^prod — ]_ _|_ -p2 ' \° V ) 

which follows from Eq. f!35p . for vanishing continuum contributions P CO nt — 0. 

Note that the coefficient ao can be obtained from the neutralino production cross 
section fT25l) . For example, for unpolarized beams, V + = V- = V = 0, 

a = aij-^=(l + 6ij). (37) 

V *ij 

Here we assume that the continuum contributions -P con t §IM to the cross section 
<7jj are already subtracted, e.g, through an extrapolation of around the reso- 
nances [22], and/or by neutralino cross section measurements at the International 
Linear Collider (ILC) [44,45]. 
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3.2 Asymmetries of the lepton energy distribution 



The longitudinal neutralino polarization is also sensitive to the Higgs interference in 
the production yU + /i~ — > X°X°. The neutralino polarization can be analyzed by the 

subsequent decays x°j ~* "li, with I — e,/i, and x°j ~~ " T ^12- ^ n Section [2751 we 
have shown that the slope of the lepton energy distribution, see Fig. [3], is proportional 
to the averaged longitudinal neutralino polarization S^ cs /P. The polarization can 
be determined by the the energy distribution asymmetry [26] 

Arr n 1 V 3 

1 „ (I + V + V-)b + (V + + V-)^ 

= nW* p 75 > ( 38 ) 

^ r res 1~ r cont 

with Acr™± = ct"± (E e > Eg) — ct™± [Eg < Eg), and n = R, L for £ = e, /i, and n — 1, 2 
for I = t. Here we have used the explicit formula for the energy distribution of 
the decay lepton £ ± ( 13T1) . with E P = E^ es ( 1331) . The continuum contributions to 
the neutralino polarization E^ ont = vanish due to the Majorana properties of the 
neutralinos, see Section [231 

The average neutralino polarization is thus solely due to the Higgs exchange, and 
receives CP-even and CP-odd contributions, proportional to b\ and b , respectively. 
In order to separate these coefficients we define the polarization asymmetries 

lin _ Aa? ± (V) - Aa? ± (-V) _ Vb, 

A± " oZW+aM-V) ~ W (1+P>o + Pco nt ' ( j 

lin Av^V) + A^j-P) _ 1 (1 + P 2 )6q 

l± " a^(V)+a^(-V) ~2 W (1 + P> + P cont ' ^ 

for equal muon beam polarizations V + = P_ = V . 

The slepton from the neutralino decay, x° — > ^ ± ^, subsequently decays into 
a neutralino and a secondary lepton. The primary and secondary leptons have 
to be distinguished from each other, for example, by using their different energy 
distributions, see Fig. [3j However, the largest part of that irreducible background 
from the secondary lepton cancels in forming the charge conjugated asymmetries [26] 

/lP ol > n \l \V°\ri _ yipol,nN _ n ^1 (at\ 

A ~ 2 ( ^ + ^" j " % (l+P 2 )a + Pco n t' ( j 

j/pol,n _ 1/ Alpol,n _ Alpo\,n\ _ 1 n U I" ^ )^0 /^x 

A " 2 ( ^ + j -2 % (l + P> + P cont - {A2) 

Due to the pure left or right coupling structure of the neutralinos to the selectrons 
and smuons, the asymmetries for the decay into e^, £ir and e^, fiL, respectively, have 
opposite sign: 

A P ol,R _ _ yipol,L A>po\,R _ _ y|/pol,L p _ f 4 0\ 
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which follows from rj^ = +1 and rj^ = — 1. The asymmetries for the decay into f% 
and T2 are always smaller than those for the decays into a selectron or smuon, 

A p T ol ' n = ^Af' R , AT l ' n = V?A'l ol ' R , 71 = 1,2, t = e,n, (44) 

with \rff\ < 1, due to mixing in the stau sector, see Eq. ( 1321) . 

The CP-even asymmetry A pol,n is due to the correlation between the longitudinal 
polarizations of the initial muons and final neutralinos, see Appendix [B] Large values 
of the CP-even asymmetry A^° ' n can be obtained when both Higgs resonances are 
nearly degenerate, and if their amplitudes are of the same magnitude. However, a 
scalar-pseudoscalar mixing in the presence of CP phases will in general increase the 
mass splitting of the Higgs bosons, and the reduced overlap of the Higgs resonances 
also reduces the CP-even asymmetry AJ ol ' n . 

The CP-odd asymmetry A'[ ol ' n vanishes for CP-conserving Higgs couplings. Sim- 
ilarly to the CP-odd polarization asymmetry .Ap°o d fl35|) for neutralino production, 

the decay asymmetry Af ol ' n is approximately maximal if the Higgs mixing is res- 
onantly enhanced. As pointed out earlier, this can happen naturally in the Higgs 
decoupling limit. 

4 Numerical results 

We analyze numerically the CP-odd asymmetry AE° od fl3"4"l) for neutralino production 
H + jjT — > XiX.2- F° r the subsequent decays, xl e &R an d x\ — * tti, we study the 
CP-even and CP-odd polarization asymmetries ^.p 01 ^, ^.p 01,1 (jH]), and Af ol ' R (|42jl . 
respectively. The feasibility of measuring the asymmetries also depends on the neu- 
tralino production cross section and decay branching ratios, which we discuss in 
detail. 

We induce CP violation in the Higgs sector by a non-vanishing phase <pA of 
the common trilinear scalar coupling parameter A t = Ab = A T = \A\exp(i(pA) for 
the third generation fermions. This assignment is also compatible with the bounds 
on CP-violating phases from experiments on electric dipole moments (EDMs) [46- 
49]. We assume a CP-conserving gaugino sector, i.e., we keep the gaugino mass 
parameters Mi, M2, and the Higgs mass parameter [i real. For the calculation of the 
Higgs masses, widths and couplings, we use the program FeynHiggs 2.5.1 [16, 17], 
see also Appendix For the branching ratios and decay width of the neutralino, 
we include the two-body decays [42] 

X°2 -> t + L + Xi + Z, Xx + Hu (45) 

with n = R, L for £ = e, /1, and n — 1, 2 for I — r. We neglect three-body decays. In 
order to enable the leptonic neutralino decays X2 ~~ > ^n, we need light sleptons. We 
parametrize their masses by ttiq and M 2 , which enter in the approximate solutions to 
the renormalization group equations, see Appendix [F] We parametrize the diagonal 
entries of the squark mass matrices by the common SUSY scale parameter Msusy = 
Mq 3 = Mjy 3 = Mf) 3 . Finally, in order to reduce the number of parameters, we 
assume the GUT relation for the gaugino mass parameters M\ = 5/3 M2 tan 2 8w- 
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Table 1: SUSY parameters for the benchmark scenario CPx- The slepton masses 
are parametrized by mo, the squark masses by Msusy- 



M H ± = 500 GeV 


tan/5 = 10 




.4 


= 1 TeV 


(f)A = 0.27T 


Msusy = 500 GeV 


[x = 400 GeV 


M 2 


= 300 GeV 


m = 100 GeV 



Table 2: SUSY masses, widths, and branching ratios for the benchmark sce- 
nario CPx, evaluated with FeynHiggs 2.5.1 [16, 17]. 



M Hl = 126.0 GeV 


m v o = 147 GeV 


m gR = 180 GeV 


B R(x° 2 - e + ~e R ) = 6.4% 


M H2 = 492.8 GeV 


m v o = 275 GeV 


m SL = 289 GeV 


BR(x° r+ff) = 23% 


M H3 = 493.6 GeV 


m x o = 405 GeV 


m fl = 178 GeV 


BR(x° - X\Z) = 9.1% 


Y H2 = 0.97 GeV 


m x ± = 274 GeV 


m f2 = 290 GeV 


BR(x§ ^ X\H X ) = 18% 


IV 3 = 0.93 GeV 


m + = 433 GeV 


m s = 277 GeV 


v l = -0.53 



We center our numerical discussion around scenario CPx, defined in Table HJ 
Inspired by the benchmark scenario CPX [50] for studying enhanced CP-violating 
Higgs-mixing phenomena, we set \A\ = 2Msusy = 1 TeV, M 3 = 800 GeV, and 
a non-vanishing phase <pA = 0.2tt. We thus obtain large contributions from the 
trilinear coupling parameter A of the third generation to the Higgs sector, both CP- 
conserving and CP-violating. In contrast to the CPX scenario, we do not need large 
values of \fi\ to obtain large CP violating effects, see discussion in Section |4~2"1 We 
choose /i = 400 GeV and M 2 = 300 GeV of similar size to enhance the branching 
ratios of the Higgs bosons into neutralinos, which are large only for mixed neutra- 
linos. We fix tan/5 = 10, since the Higgs boson decays into neutralinos are most 
relevant for intermediate values of tan (3. Smaller values of tan j3 favor the tt decay 
channel, while larger values enhance decays into bb and tt. We give the masses of 
the Higgs bosons, the charginos, neutralinos, light sleptons, and the widths of the 
Higgs bosons for scenario CP% in Table [2} where we also list the branching ratios 
for the decaying neutralino. We choose longitudinal muon beam polarizations of 
V + = T-L = V = ±0.3, which should be feasible at a muon collider [37]. 

4.1 yfs dependence 

For the scenario CPx, we analyze the dependence of the asymmetries and the cross 
section on the center-of-mass energy y/s. The CP-even and CP-odd observables 
exhibit a characteristic y/s dependence, mainly given by the product of Higgs boson 
propagators Am), see Eq. ( 1231) . A muon collider will have a very precise beam 
energy resolution, and thus enables detailed line-shape scans. 

In Fig. H](a), we show the CP-odd polarization asymmetry ^4p°o d for neutralino 
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production — > X1X2 as a function of y/s around the heavy Higgs resonances 
if 2 and if 3. At the peak value, y/s = (M H2 + M Hs )/2 ps 493 GeV, the interference 
of the two nearly degenerate Higgs bosons is maximal, leading to an asymmetry 
of up to Apf od = 30%. The asymmetry measures the difference of the neutra- 
lino production cross section ai2(V) for equal positive and negative muon beam 
polarizations V = ±0.3, which we show in Fig. S(b). We also observe that the 
splitting of the two resonances is increased in the presence of CP-violating phases. 
For <pA = 0.27T, the two resonances are clearly visible in the line shape of the cross 
section a 12, whereas it assumes a single resonance form for 0^ = 0, where the Higgs 
bosons are extremely degenerate, see Fig. H(b). 

The Higgs boson interference in neutralino production also leads to CP-even and 
CP-odd contributions to the average longitudinal neutralino polarizations. In order 
to analyze the X2 polarization, we discuss the CP-even asymmetry A^ ol,R and the 
CP-odd asymmetry Af ' of the leptonic energy distributions for the neutralino 
decay. For simplicity, we discuss only the decay into a right selectron X2 ~^ e &R- 
The same asymmetries are obtained for the decay into a right smuon X2 ~~ > Wr- 
For the neutralino decay into a tau, X2 ~ * T ^i> t ne corresponding asymmetries are 
obtained from Eq. (J44J). 

For the decay X2 ~~ * e ^Ri we show the y/s dependence of the CP- and CPT-even 
asymmetry ^p 01 ^ in Fig. [1(c) for <\> A = and <p A = 0.2%. The CP- and CPT-odd 
asymmetry Af ol,R is shown in Fig. H](d). The phase 4>a tends to increase the mass 
splitting of the Higgs resonances. Their overlap is now reduced, leading in general 
to a suppression of the CP-even asymmetry A^ ol ' R , in particular at the mean energy 
of the resonances y/s = (Mh 2 + Mh 3 )/2, see Fig. 11(c). On the contrary, the larger 
Higgs splitting increases the CP-odd asymmetries Af ol,R and -4.p°o d - 

All asymmetries for production and decay vanish asymptotically far from the 
resonance region. The continuum contributions from selectron and Z exchange 
to the difference of the cross sections and to the average neutralino polarization 
cancel in the numerator, but contribute in the denominator of the corresponding 
asymmetries, see their definitions in Sectional 

In the following Sections, we analyze the dependence of the production cross 
section and the asymmetries on |A| and 4>a, and finally on M2 and /x, fixing all 
remaining parameters to those of scenario CP%- We fix the center-of-mass energy 
to y/s = (Mh 2 + Mff 3 )/2, where we expect the largest CP-odd asymmetries Apf od 
and A'p o1 ' r , see Figs. H^a) and (d), respectively. For consistency, we also choose 
y/s = (Mh 2 +Mh 3 )/2 for the discussion of the CP-even asymmetry ^.p 01 ^, although 
it is generally suppressed at this value if CP is violated. 



4.2 \A\ and (J)a dependence 

We analyze the dependence of the CP-asymmetries on the phase 4>a of the trilinear 
coupling A, which is the only source of CP violation in our study. The CP-odd 
asymmetries, AZ°l d and A^ ol,R , see Fig. |5(a), are approximately maximal, if the 
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Figure 4: ^fs dependence of (a) the CP-odd production asymmetry «4.p° od , Eq. ([Ml), and 
(b) the cross section a\2 for neutralino production /i + pT — > x±X2- For the subsequent 
decay, x° ~* e &Ri lu ( c ) the CP-even polarization asymmetry AT l ' R , Eq. (gl]), and m 
(d) the CP-odd polarization asymmetry AT° l,R , Eq. (T42l) . for the SUSY parameters as 
given in Table [TJ The longitudinal beam polarizations are V+ = V- =V = +0.3 (solid), 
in (b) V = -0.3 (dashed), and in (b,c) V = +0.3 with <f> A = (dot-dashed). 



mixing of the Higgs states is resonantly enhanced. This is naturally achieved when 
the diagonal elements m 2 H — S##(s) and m 2 A — T<aa(s) of the Higgs mass matrix 
M, Eq. ([H]), are equal, provided the corresponding imaginary part is small. In 
our scenario, where decays into heavy squarks are not kinematically allowed, this 
condition is roughly fulfilled for <pA — 0.27T. We interpret this condition as a level 
crossing of the CP-eigenstates H and A, when m 2 H — m 2 A — Re[£HH"(s) — Saa(s)] 
changes sign [10]. The mass difference of the physical Higgs boson masses, however, 
is typically increased by the H-A mixing, as can be observed from Fig. 0(c). A 
splitting of the order of the Higgs widths Th 23 , shown in Fig. E](d), leads to large 
absorptive phases, which are necessary for the presence of CPT-odd observables. The 
increased Higgs mass splitting for non-vanishing phases leads, however, in general 
to lower peak cross sections auiV), which we show in Fig. Mb), both for positive 
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Figure 5: Phase dependence of (a) the CP-odd polarization asymmetry «4.p° od (solid), 
Eq. (|34"|) . for neutralino production fi + [i~ — > X1X2 a * v 7 * = (Mh 3 + Mff 3 )/2, and for 
the subsequent decay X2 ~~ * the CP-even polarization asymmetry Ae° l,R (dashed), 
Eq. (I4ip . and the CP-odd polarization asymmetry 

X pol ' R (dot-dashed), Eq. (02]). In 
(b) the cross section a±2 for neutralino production /U + /i~ — > X1X2 with longitudinal beam 
polarizations V+ = V- = V = +0.3 (solid), and V = —0.3 (dotted). In (c) the Higgs 
masses Mn i and (d) the Higgs widths r# i; for i = 2 (solid), and i = 3 (dashed). The 
SUSY parameters are given in Table [TJ 



and negative beam polarizations V = ±0.3. 

The asymmetries and cross section for negative <pA can be obtained from symme- 
try considerations. Since the complex trilinear coupling A is the only source of CP 
violation in our analysis, the CP-odd asymmetries -4.p°o d and ^4' e pol ' R must be odd 
with respect to the transformation <f) A — > —rf'A, while the CP-even asymmetry A^ ol ' R 
must be even. Consequently, the cross section transforms as a^iV) — > cr 12 (— V). 

In Fig. El we show contour lines of the cross section and the asymmetries in 
the 4>a~\ A I plane. The largest CP-odd asymmetries A^°l d and A'^ ol ' R are obtained 
for \A\ ~ 2Msusy = 1 TeV. For larger values of \A\, the lighter stops become 
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kinematically accessible and H 2 decays dominantly into tfti pairs, which leads to 
a suppression of the neutralino production cross section. We therefore restrict our 
discussion to \A\ < 1.2 TeV. 

As we have observed in Fig. 0(a) of the preceding paragraph, the CP-even asym- 
metry A p e ol ' R , Fig.EKc), is in general larger in the CP-conserving limit. The maxi- 
mum of the asymmetry is also obtained for A « ±800 GeV. However, this is rather 
coincidental, and is due to the exact degeneracy of the Higgs bosons H and A. 

Note that large resonant mixing is possible without requiring large values of \fi\. 
The CP- violating scalar-pseudoscalar self energy transitions are proportional to the 
amount of CP violation in the squark sector, described by the quantity 

3 Im ^) (46) 
167r 2 m 2 - — m 2 ~ 

h h 

with f = t,b [6,8,28]. However we obtain large H-A mixing for moderate values of 
fi, since, in the Higgs decoupling limit IhiShh(s) ~ IulSaa(s) for energies below the 
threshold of heavy squark pair production. Therefore, the conditions for maximally 
resonance enhanced mixing discussed in this section may be fulfilled for moderate 
values of T*ha(s). 



4.3 fi and Mi dependence 

The couplings of the Higgs bosons to the neutralinos strongly depend on the gaugino- 
higgsino composition of the neutralinos, which are mainly determined by the values 
of \x and M 2 . For neutralino production —> XiX^ we show the CP-odd po- 

larization asymmetry .4p° od flMl) in the /i-M 2 plane in Fig. 0(a). For V = ±0.3, 
the maximum absolute value of the asymmetry would be A 1 ^™ 3 * ~ 55%, as follows 
from Eq. (J36l) . We observe in Fig. [71(a) that the asymmetry reaches 40% near the 
neutralino production threshold, where the coefficient a\, see Eq. (1181) . receives large 
spin-flip contributions. For smaller \i and M 2 the Higgs boson widths are increased, 
since decay channels into light neutralinos and charginos open. This results in a 
larger overlap of the Higgs resonances, which reduces the absorptive phases, and 
consequently suppresses the CPT-odd asymmetry «4£° od . In the upper left corner of 
Fig. [7(a), the asymmetry changes sign due to a level crossing of the two neutralinos 

xl and xl- 

In Fig. E(a), we show the cross section <7i 2 for neutralino production — > 
XiX2- m the mixed region ~ M 2 , where the Higgs-neutralino couplings are larger, 
the cross section reaches up to cr 12 w 1500 fb. In addition, since H 2 and H s are mixed 
CP-eigenstates, there is no p-wave suppression. However, due to the Majorana 
nature of the neutralinos, the continuum contribution from fx and Z exchange to 
the cross section is p-wave suppressed [51]. It is thus negligible near threshold, and 
reaches 150 fb only for fi < 150 GeV. In Fig. E](b), we show the branching ratio for 
the neutralino decay x 2 ~* e+ &R- The decay fraction is reduced in the upper right 
corner, since the channels x 2 ~~ ¥ xH^i an d X 2 ~^ Xi-Hi, open. In particular, the 
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(c) 4>a/^ (d) ^a/tt 

Figure 6: Contour lines in the <^4-|.A| plane for (a) the cross section for neutralino 
production — > X1X2) O 3 ) the CP-odd polarization asymmetry -4.p°Q d , and for the 
subsequent decay X2 ~^ e ^R> ( c ) the CP-even polarization asymmetry ^l| o1 '^, Eq. (|4ip . 
and (d) the CP-odd polarization asymmetry X P ° 1,R , Eq. (g2]), at ^/s = (M H2 + M H[1 )/2 
with longitudinally polarized beams "P_ = "P+ = ±0.3. The SUSY parameters are given 
in Table H 



branching ratio into the lightest Higgs boson can be BR(x2 —> Xi^i) > 60%, for 
M 2 > 400 GeV. 

For the neutralino decay X2 ~^ e ^Ri we show the polarization asymmetries A^ ol ' R 
and Af ol ' R in Fig. [9](a) and (b), respectively. As discussed before, the CP-even 
asymmetry A^ ol,R is suppressed by CP-violating effects due to the smaller overlap 
of the resonances at a/s = (M# 2 + M# 3 )/2. Therefore we only find large values 
of A^ ol,R for light neutralinos and charginos in the lower left corner of Fig. 0(a), 
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Figure 7: Neutralino production ji + fi — > X1X2 at = (Mff a + Mh 3 )/2 with longitu- 
dinally polarized beams T 3 - = = —0.3. Contour lines in the [i—M-i plane for (a) the 
CP-odd production asymmetry ^4p°o d AMD , and (b) the significance S^ oA (|H.70p . with 
C = 1 fb _1 , for the SUSY parameters as given in Table [TJ The corresponding cross 
section a(fi + fi~ — > X1X2) 1S given in Fig. [HJ The shaded area is excluded by requiring 
m + > 103 GeV. 
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Figure 8: Contour lines in the \1-M1 of (a) the cross section o\2 = a(fi + ^ — ► X1X2) 
at \/s = (Mh 2 + M# 3 )/2 with longitudinally polarized beams V- = V+ = 0.3, (b) the 
branching ratio BR(%2 —* e+ e^), for the SUSY parameters as given in Table [TJ The 
shaded area is excluded by requiring m ± > 103 GeV. 



where the larger Higgs widths counter the effect of the larger Higgs mass difference. 
On the contrary, in that region the CP-odd asymmetry A!£ oX,R is reduced due to 
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Figure 9: Neutralino production fi + fi — > X1X2 an d decay X2 ~^ e &R a t y/s = (Mff 2 + 
Mjj 3 )/2 with longitudinally polarized beams V~ = = ±0.3. Contour lines in the 
\1-M2 plane for (a) the CP-even polarization asymmetry Ae° l,R , Eq. (|4ip . (b) the CP- 
odd polarization asymmetry A'i° l ' R , Eq. flM]), and (c) the significance Sf ol ' R , Eq. P02|) . 
and (d) the significance 5e pol ' R , Eq. (|H.73|) . with the effective luminosity £ c ff = 1 fb -1 . 
The SUSY parameters are given in Table [TJ The neutralino production cross section, 
a(fi + fi~ — > X1X2) an d branching ratio BR(x!] — > e+ ^~R) are shown in Fig. El The shaded 
area is excluded by requiring m ± > 103 GeV. 

Xi 



smaller absorptive phases. Finally, at threshold the longitudinal polarization of the 
neutralino Ej? es (|T5l) vanishes, and thus also both decay asymmetries, as follows from 
Eqs. (DSD and flU}. 

The significance of the CP-odd polarization asymmetry, defined in Appendix [Hi 
reaches S'? 01,11 ~ 2, see Fig. M^c), and thus the measurement of this asymmetry will 
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be challenging. Nonetheless, only the asymmetry Af ol ' n allows to measure the CP- 
odd contribution b to the longitudinal neutralino polarization Sj? es ( 1151) . However, 
taking the other leptonic neutralino decay modes into account, in particular ~^ 
rfi, to analyze A!? 01 ' 1 , results in larger significances^]. 

4.3.1 Neutralino decay into a stau-tau pair 

The CP-even and CP-odd neutralino polarization asymmetries can also be measured 
for the neutralino decay into a tau, X2 ~~ * r ^i- Due to mixing in the stau sector, 
the asymmetries for the decay into a tau are generally smaller than those for the 
decay into an electron (or muon), *4.p o1,1 = rfrA^ 1 ' 1 and A'f o1,1 = r]lAf ol ' R , with 
\t)t\ ^ 1) see Eq. 044p . In Fig. fTDT b). we show the contour lines of the reduction 
factor 7]\ (1321 . In the following, we discuss the CP-even asymmetry A^.° 1,1 only. 
A similar discussion holds however qualitatively also for the CP-odd asymmetry 
A'f Note that a measurement of the r asymmetries is more involved due to r- 
reconstruction efficiencies, which we however neglect in the following for simplicity. 

The CP-even asymmetry A^.° 1,1 is shown in Fig. flDT c). The significance for 
measuring the asymmetry also depends on the cross section for production and decay, 
and thus on the different leptonic branching ratios of the neutralino. The neutralino 
decay into a tau dominates for M 2 < 200 GeV, see Fig. HOT c). whereas the branching 
ratio into an electron can attain more than BR(x2 ~^ e+ ^j?) = 16%, for \x < M2, see 
the contour line in Fig. [5](b) . We take account of this interplay between the size of 
the branching ratios and the asymmetries by comparing their statistical significances 
gpoi,R anc j 5poi,i^ w hich we define in Appendix |Hl Eq. (IH.72I) . The significances 
quantify the feasibility of measuring the asymmetries. Both significances can be as 
large as 5, however in different regions of the /x-M 2 plane, compare Fig. M^c) and 
Fig. ITuTd). respectively. 



3 The decay asymmetries Ae° 1,L , Aji ol ' L , and At^ p ° 1 ' 2 are only accessible for [i > 500 GeV 
and M2 > 200 GeV in our scenario. They are not relevant for our discussion, since the correspond- 
ing branching ratios of \\ are only as large as a few percent. 
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Figure 10: Neutralino production 

\fs = (Mh 2 + Mh 3 )/2 with longitudinally polarized beams V- 

T + f- 



>~0 



tt\ 1 at 

V+ = ±0.3. Contour 

lines in the /j— M2 plane for (a) the branching ratio BR(x2 — ► t+ ^i)j (b) the factor 77* (|32p . 
(c) the CP-even polarization asymmetry ^l? 01 ' 1 , Eq. (|4ip . and (d) the significance S^ ' 1 , 
Eq. (|H.72h . with the effective luminosity £ e g = 1 fb" 1 , for the SUSY parameters as given 
in Table [TJ The shaded area is excluded by requiring rn x ± > 103 GeV. The cross section 

a(fi + fi~ — > X1X2) is shown in Fig. El 
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5 Summary and conclusions 



We have analyzed neutralino production and their leptonic decays at the muon col- 
lider with longitudinally polarized beams. We have defined polarization asymmetries 
to study the interference of the heavy neutral MSSM Higgs bosons with CP viola- 
tion, radiatively induced by the common phase 4>a of the trilinear scalar coupling 
parameter. For nearly degenerate neutral Higgs bosons, as in the Higgs decoupling 
limit, the CP violating Higgs mixing can be resonantly enhanced, which allows for 
large CP violating effects. 

For neutralino production, we have defined a CP-odd asymmetry of the cross 
section for equal positive and negative muon beam polarizations. The CP-odd pro- 
duction asymmetry is sensitive to the CP-phases in the Higgs sector, and also re- 
ceives large contributions from absorptive phases of the Higgs propagators. In a 
numerical study, we have obtained large values of the production asymmetry up to 
40% for equal beam polarizations of V = 0.3. For neutralinos with mixed gaugino- 
higgsino character, the production cross section can be as large as 1500 fb. Thus 
the asymmetry can be measured with high statistical significance. 

We have shown that also the neutralino polarization depends sensitively on the 
Higgs interference. For the subsequent leptonic decays of the neutralino, we have 
analyzed two asymmetries of the energy distributions of the final leptons e, \x and 
t. The decay asymmetries probe the CP-odd and the CP-even contributions to 
the longitudinal neutralino polarization, respectively. The decay asymmetries are 
complementary to the production asymmetry, since they strongly depend on spin- 
correlations. The CP-even asymmetry is due to a correlation between the longitudi- 
nal polarizations of the initial muons and the final neutralinos. Being CP-even, the 
asymmetry reaches 25% for vanishing CP-phases, and is reduced in the presence of 
CP-phases. The CP-odd asymmetry is due to the spin correlations in the neutralino 
production and decay process. Similarly to the CP-odd asymmetry from the pro- 
duction, this decay asymmetry is approximately maximal if the scalar-pseudoscalar 
Higgs mixing is resonantly enhanced, which appears naturally in the Higgs decou- 
pling limit. The decay asymmetries yield additional information on the CP nature 
of the Higgs resonances, and complement the production asymmetry. The asym- 
metries thus allow a systematic study of the interference and mixing effects of CP 
violating neutral Higgs sector at the muon collider. 
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Appendix 



A Note on the effective Higgs couplings 

In Section 12.11 we have denned the effective couplings of the Higgs bosons to muons 
and neutralinos. They are obtained by rotating the tree level Higgs couplings by the 
matrix C. This matrix includes the leading radiative self-energy corrections, and is 
defined by diagonalizing the Higgs propagator matrix and the Higgs mass matrix 

A D (p 2 ) = CA(p 2 )C~\ M D (p 2 ) = CM(p 2 )C~\ (A.l) 

respectively, at fixed momentum squared p 2 = Mfj 2 . The weak momentum depen- 
dence of the mass matrix M is neglected, as in the Weisskopf-Wigner approxima- 
tion [39]. Due to the absorptive parts of the transition amplitude, the matrix C is in 
general non-unitary. Thus the transformed basis of the approximate^ Higgs boson 
fields {Hi, H 2 , H 3 }, is non-orthonormal [10]. As a consequence, there exists a dual 
basis {Hi, H 2 , H 3 } obtained by the matrix C = C~ l . The corresponding states 
satisfy the orthogonality relations (Hk\Hi) = 5m- This leads to different transfor- 
mations of the tree-level Higgs couplings to initial and final state fermions with C 
and C, relations (jHJ) and (jBJ), respectively. This follows, since the amplitude ( JTUj) for 
neutralino production can also be written in the general form 

= t^c^cac^ct^ 

= tMc^AdCTM = T^A d T^, (A.2) 

where r&> and are the one-particle irreducible Higgs vertices to muons and neu- 
tralinos, respectively. Eq. ( 1A.2I) defines the effective one-particle irreducible Higgs 
vertices for initial and final fermion states, 

rffi = QjFf, (A.3) 

r& = ciwrw (A.4) 

which transform with C and C, respectively. If the phases of the Higgs boson states 
are chosen appropriately, the matrix C can be made complex orthogonal [52], which 
implies C = C . 

B Correlation between initial and final longitudi- 
nal polarizations 

In this appendix, we analyze the correlation between initial and final longitudinal 
polarizations in neutralino pair production in annihilation via Higgs boson 

exchange in the simplified case of a CP conserving Higgs sector. 

4 We call the Higgs basis {Hi, H2, H3} approximate, since it corresponds to Md(p 2 ) with fixed 
momentum squared p 2 = Mfj , assuming a weak momentum dependence in the resonance region. 



24 



Helicity and CP eigenstates 

Out of four possible spin/helicity states of a fermion pair only those with J z = in 
the center-of-mass system (CMS) interact with the Higgs bosons. Here J denotes 
the total angular momentum and z the direction of the momenta of the fermions. 
Since in the CMS the orbital angular momentum of the fermions is orthogonal to 
their momenta, L z = 0, their total spin S satisfies S z = J z — L z = 0. The fermion 
interacting states are thus \LL)ff> and \RR)ff, where L and R denote the helicities 
of the fermions //' = fi + n~ , XiXj- These states are linear combinations of states 
IS", S z ) ffi with spin S — and S = 1, 

1 

3' 



\LL) ff/ = —(\l,0)-\0,0)) ff , (B.5) 



\RR)ff = ^(M + M)ff- (B.6) 

For a fermion-antifermion system their spin S is related to their CP quantum number 
by CP = rjff(— 1) 5+1 , with = 1 for a Dirac fermion-antifermion pair and tj x o x o = 

= e 2 " 1 " 13 for a pair of neutralinos. The relative CP phase factor rjij is real in our 
analysis since the neutralino sector is CP-conserving, with = 0, tt/2. The CP- 
even and CP-odd muon states 

|CP+) MV -= |1,0) M+At -, (B.7) 
|CP-) M+M - =i|0,0) M+M -, (B.8) 

and neutralino states 

|CP+) x o x o = ( coscrjj|l, 0) + % sin cry 1 0, 0)) x o x o, (B.9) 
|CP— ) x o x o = (— sin <7jj|l, 0) + i coscrjj|0, 0)) x o x q, (B.10) 

can be expressed as a linear combination of helicity states inverting Eqs. (IB. 51) and 
( IB.6p . Analogously, the helicity states are linear combinations of the CP-even and 



CP-odd states. For the muon-antimuon pairs we obtain 

1 

1 

7! 



\LL)^- = ^(|CP+> + i|CP-» M+M -, (B.ll) 
\RR)^ +IM - = 4s(l CP +) " ^l cp -))^-- (B.12) 



Transition amplitudes 

Assuming a CP conserving Higgs sector implies that, in our Higgs mediated neu- 
tralino production process, a CP-even state ICP-I-)^^- can only produce 
CP-even Higgs bosons, which in turn decay into the CP-even neutralino state ( IB. 91) 
with a real amplitude a. Analogously, a CP-odd state leads to a CP-odd 
neutralino state (IB.lOp with an amplitude (3. 

An initial state with right handed polarized muon and antimuons (IB. 121) . will 
interact to produce the neutralino state 

IflflW- -> v / 2A^[a|CP+) -z/5|CP-)] x o x o, 

= N' R [e 1 ^ (a + f3)\RR) + {a - (3)\LL)] x o x o, (B.13) 
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where M' R is a normalization factor. Here we have used the explicit form of the 
neutralino states (IB. 91) and (IB.lOp and have inverted Eqs. (1B.5P and (IB. 61) . 
Similarly, 

\LL)n+n- ^Mi[e ia »{a- P)\RR) + e^(a + P)\LL)] x o x o. (B.14) 

From Eqs. ( 1B.13I) and (IB. 141) follows that if either the CP-even or CP-odd ampli- 
tudes, a or j3, respectively, vanish, then so does the neutralino polarization, since in 
this case the absolute value of the coefficients on the r.h.s. of Eqs. (IB. 13|) and (IB. 14[) 



are equal. Note that the neutralino polarization depends on the relative signs of the 
CP-even and CP-odd amplitudes since it arises from their interference. 

From Eqs. (IB. 13[) and (IB. 14[) also follows that if the muon beams are not longi- 
tudinally polarized, which implies that the initial state has equal proportions of left 
and right handed states ( IB. Ill) and (IB. 121) . the final neutralinos will also be 

unpolarized. 

Concluding, in the CP conserving Higgs sector the neutralino polarization is 
sensitive to the relative sign of the transition amplitudes, and thus to the product 
of couplings, and can only be non vanishing if the muon beams are longitudinally 
polarized. This implies that the correlation between initial and final state longitudi- 
nal polarizations depends on the interference of transition amplitudes mediated by 
Higgs bosons of different CP parities. 

In the more general CP violating Higgs sector studied in this paper, the above- 
mentioned correlation between initial and final polarizations leads to the CP-even 
asymmetry 

A Po\,n gjjj In additioil) CP-odd (and CPT-odd) effects lead to the 
CP-odd asymmetry A!f ol ' n ( 1421) . However, since this asymmetry is not due to a 
correlation between initial and final state polarizations, it can be non-zero even for 
vanishing beam polarizations. 



C Lagrangians for non-resonant neutralino pro- 
duction and leptonic decay 

The non-resonant neutralino production (CQ) proceeds via Z° boson exchange in 
the s-channel, and smuon fiL.R exchange in the t- and w-channels, see the Feynman 
diagrams in Fig. [21 The interaction Lagrangians for neutralino production and those 
for its leptonic decay Xi ~^ ^l,r, with £ = e,/i are [41,53] 

^°x°x° = \z^ Z ltnO^P L + 0:?P R }x% *,j = l,...,4, (C.15) 

J Zi Cub u\y 

&z°w = 9 —Z° u ^[L^P L + R^P R ]fi, (C.16) 

COS U\y 

&t®> = grt j £PR$h + 9fg£PLX j £R + h.c.. (C.17) 
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with Pl,r = (1 : F75)/2. In the photino, zino, higgsino basis (7, Z, H®, H®), the 
couplings are [41,53] 

-\ [(N i3 N* 3 - N i4 N* 4 ) cos 2(3 + (N i3 N* 4 + N t ,N* 3 ) sin 2(3] 



O L 



O 



'R 



1.1 



-°4 



1 

2 



sin 2 W 
1 

cos 9w 
i sin 



sin 2 # 



(T 3£ - e e sin 2 6 w )N j2 + e e sin 6 w N jl 



tan6 w N% - N* x 



(C.18) 
(C.19) 
(C.20) 

(C.21) 

(C.22) 



with ee and T 3 £ the electric charge and third component of the weak isospin of the 
lepton £, the weak mixing angle 8w, the weak coupling constant g = e/ sin 6w, 
e > 0, and the ratio tan/5 = v 2 /vi of the vacuum expectation values of the two 



neutral Higgs fields. The neutralino couplings to the Z boson, O i ■ ' , and to the 



smuons, / ■ , contain the complex mixing elements JVj,-, which diagonalize the neu- 



1.1 



tralino matrix N* a Y a pNl k = m Xi 5ik [1], with neutralino masses m Xi > 0. 

Mixing can safely be neglected for the scalar leptons of the first two generations, 
£ — e,jl. For the neutralino decay into staus x® — > f n r, n 
mixing into account, see Appendix O 



1,2, we take stau 



D Density matrix formalism 



We use the spin density matrix formalism of [40,41] for the calculation of the squared 
amplitudes for neutralino production ([T]) and decay The amplitude for neutra- 
lino production via resonant Higgs exchange, Eq. flD.23p . depends on the helicities 



A-t of the muons /j? and the helicities A,;, Xj of the neutralinos x° 



T 



A(H k ) v(p, + ,\ + )[c»^P L 



c H R ^P R )u{p,-A~) 



fi(PxS» A i) [*L XlXip L + 4 kX ^P R ) v(p x o, \)} . (D.23) 



We include the longitudinal beam polarizations of the muon-beams, V- and V. 
with — 1 < V± < +1 in their density matrices 

1 



Ptx> 



<$A_A'_ 



(D.24) 



(D.25) 



where r 3 is the third Pauli matrix. The unnormalized spin density matrix of Xi Xj 
production and x° decay are given by, respectively, 



p 

Px,X'. 



E+ - rpP rpP* 

Px+X' Pa_ X'_ 1 AiAj :A+A_ 1 X t X' ' X' + X'_ 



[D.26) 



Xi,X+,X',,X-X'_ 



27 



p°, x . = T%T°. (D.27) 
The amplitude squared for production and decay is then 

\T\ 2 = |A(x^E^4v (D.28) 

with the neutralino propagator 

A(X°) = ~ =-. (D.29) 

The spin density matrices ( ID. 261) and ( ID. 271) can be expanded in terms of the Pauli 
matrices T a 

3 

P P x jK = V; f + E^Sp, (D.30) 

a=l 



^ = fe^ + E^^ (D.31) 

a=l 

where we have defined a set of neutralino spin vectors s° - In the center-of-mass 
system, they are 

S y = (0; 1,0,0), sV = (0; 0,1,0), A" = — (|0 x o|; 0, 0, £ x o). (D.32) 

We have chosen a coordinate frame such that the momentum of the neutralino Xj 
is given by 

^o = (^ x o;0,0,b- x o|), (D.33) 

with 

s + m 2 — m 2 f\ 



27i ' 2^ 



E x o = ^j, ^, l^o | = (D.34) 



= X(s, m 2 , m 2 ), (D.35) 



and the triangle function 

with A (a;, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz). 

Inserting the density matrices (1D.3Q|) and (ID.31[) into (1D.28[) . gives then the 



amplitude squared in the form of Eq. 012|) . 

E Continuum amplitudes and contributions 

The amplitudes for non-resonant Z and P>l,r exchange are 

T[ iXjX+x _(s,Z) = ^-A s (Z)[v(p ll+ ,X + )r(L ll P L + R ll P R )u(p^,X- 

COS " 
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X 



T 



A j A j A_|_ A. 



t p 



.(t,/J.L) 



x 



utp&^MO'tPL + Ol R P R )v(p x o, X 



-g 2 f*CA\ji R ) v( Pti+ ,\ + )P L v(p x °, \) 



X 



u(p x o, Xj)Pmip,j.-, A- 



pj a,-a+a_ («, Ml) 



g 2 rt*riA u (fi L ) v( PlM+ ,X + )P R v(p x o, X, 



X 



A<jAj A-|_A_ 



{u,fi R ) 



u(p x o, Ai)Pi«(p M -, A_ 
g 2 fX*f*A u (fi R ) [tJ(p M+) A + )P^(p x o, A,) 



X 



with the propagators 



A S (Z) 



s — m 



2 ' 



AW) 



£ — m~ 



A<W) 



u — m 



(E.36) 



(E.37) 



(E.38) 



(E.39) 



(E.40) 



(E.41) 



and i = (p At - — £> x o) 2 and w = (p^- — p x o) 2 . We neglect the Z-width in the propagator 
A S (Z) for energies beyond the resonance. The Feynman diagrams are shown in 
Fig. [2l For e + e~ collisions, the amplitudes are given in [41]. 

The continuum contributions P con t are those from the non-resonant Z and jj>L,R 
exchange channels. The coefficient P con t is independent of the neutralino polariza- 
tion. It can be decomposed into contributions from the different continuum channels 



P(ZZ) + P(Z/^) + PGifl&O + P(%) + P(fiLp>L), 



(E.42) 



with 



P(ZZ) 



P(^j 



P(P>rP-r) 



\^{Z)\\E^ R + L^ L )E- t 

2 m 



COS 4 

x{|Oy| 2 (P x oP x o + g 2 cos 



■peO;' R ) 2 -(/mO" fl ) 2 ]m x om x o 



(E.43) 



cos 2 9 



-R,c R EiRe\A s {Z) 



+{^& R )f*:f*o:* + a u *(^)/5/5*o:;^)(p x op x o + g 2 cos 2 ^ 

-{^ u {Jk R )f^f^ R - A u \^ R )f R f R ;o'^)2E b qcosd\ }, (E.44) 

^^{i/3H/5i 2 x 
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(\A t (fi R )\ 2 + \A u (fl R )\ 2 )(E^E x o + q 2 cos 2 9) 
-(IA'OIr)] 2 - \A u (jl R )\ 2 )2E b qco S e 

-^e{(/,t) 2 (^) 2 A"(/i R )A'*( / i /? )}2m x om x o}. (E.45) 

To obtain the quantities P(Z£il), P(£ll£I'l) one has to exchange in ( IE. 441) and ( IE. 451) 

A*(/jr) —> A t (p, L ), A u (fx R ) - A u (fl L ), c R ^c L 

Rn^> Ly,, O i:) R -> , fj% -> f^, -> (E.46) 

The longitudinal beam polarizations are included in the weighting factors 

c L = (l-V.)(l + V+), c R =(l+V-)(l- V+). (E.47) 

For e + e~ collisions, the P terms are also given in [41], however they differ by a factor 
of 2 in our notation (1121) . The continuum contributions S" ont (1T31) to the neutralino 
polarization can also be found in [41], also differing by a factor of 2. However, due 
to the Majorana character of the neutralinos, the continuum contributions S" ont 
are forward-backward antisymmetric, and vanish if integrated over the neutralino 
production angle [43], see Eq. (J3"3"]l . 



F Neutralino decay into leptons 

The expansion coefficients of the decay matrix (ID. 311) for the neutralino decay into 
right sleptons x°j ~ * w hh £ = e, fj,, are 

D = ^\m 2 (m 2 x o-m}), (FAS) 

X a D = +g 2 \fi\ 2 m x o(s xy p h ). (F.49) 

For the decay into the left sleptons Xj ~~ > i - + ^Z' ^ = e,/i, the coefficients are 

D = yl/^l'Ko-mf), (F.50) 

= -9 2 \fa\ 2 m x o(s a x o- Pil ). (F.51) 

For the decay into the stau Xj ~ > T+ Tjt\ k — 1, 2, one obtains 

D = y(l<f + l^| 2 )K«-0' (F.52) 

= -9 2 (\al/-\bl/)m x o(s x y Pei ). (F.53) 

The coefficients SJ, for the charge conjugated processes, x°j ~~ > £~£ + , is obtained by 
inverting the signs of (fR49l) . (fRoTI) and (lR53l) . 
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With these definitions we can rewrite the factor E|,, that multiplies the longitu- 
dinal neutralino polarization £p in (1121) . 

£d = Ve±^(E e -E e ), (F.54) 



where we have used 



m x ,{s\ ■ Pt ) = --^(E e - E t ). (F.55) 



In order to reduce the free MSSM parameters, we parametrize the slepton masses 
with their approximate renormalization group equations (RGE) [54] 

m | = m 2 + m^ + 0.23M 2 - m| cos 2/3sin 2 (F.56) 

m) = m 2 + m^ + 0.79M 2 2 + m|cos2/3(-^ + sin 2 9 W ), (F.57) 

ml t = ml + m ] + 0.79M 2 + ~m| cos 2/5, (F.58) 

with mo the common scalar mass parameter the GUT scale. 



G Stau-neutralino couplings 

For the neutralino decay into staus x° ~~ > TkT, we take stau mixing into account and 
write for the Lagrangian [55]: 

^rr Xi = ghHalPn + biPM + h -c- , k = 1, 2; i = 1, . . . , 4, (G.59) 

with 



« 



J, = (^L)*A T - 6«=(^D*^ (n = L,R), (G.60) 



with 7££ n given in (IG.67p . The couplings and are defined by Eqs. (10.211) and 
(IG.22p . respectively, and 



h L Tj = {h«Y = -Y T (N* 3 cos f3 + N* 4 sin/3), (G.62) 

Yr = ro ^ R , (G.63) 
V 2 ttivf cos p 

with mvi/ the mass of the W boson, m T the mass of the r-lepton and N the neu- 
tralino mixing matrix in the 7, Z, Hf, H® basis. The masses and couplings of the 
r-sleptons follow from the tl — mass matrix 



= ~(rl fS) , (G-64) 




)l 
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with m~ and m~ given by Eqs. (1F.56[) and (IF. 57ft replacing mj by m 2 ., and 



Af = A T — fi* tan /3, 
ip f = axg[Af], 



(G.65) 
(G.66) 



with A T = A the (common) trilinear scalar coupling parameter. The r mass eigen- 
states are (ti,t 2 ) = (tl,tr)71 t , with 



The mixing angle is 



COS#f 



-ra r |Af I 



e 1 ^ cos 0f sin 6^ 
— sin Of e~ % ^ f cos Of 

sin 0~ 



(G.67) 



2 2 

m~ — m| 



Jrr^\Kf | 2 + (m^ - mfj 2 



and the mass eigenvalues are 

1 



m- 



m~ + m? ) =F A / (mf - m? ) 2 + 4m 2 |A f | 2 



(G.68) 
(G.69) 



H Statistical significances 



We define the statistical significance of the asymmetry Apf od ( 1341) by 

^prod = l^prodl V^ij^, 

where £ denotes the integrated luminosity and 



a 



(H.70) 



:h.7i) 



is the mean value of the neutralino production cross section a^, (1251) . for both 
equal muon beam polarizations V and — V . The significances for the polarization 
asymmetries A v i ° [,n ( j4"Ti) and Af° ' n ( 1421) are defined by 



pol.n 



14 



pol,n 



4^BR(x 



and 



5 



/pol,n 



1-4' 



/pol,n I 



%BR(x^«-)4ff, 



(H.72) 



(H.73) 



respectively, where £ cff = e™£ is the effective integrated luminosity, with e" the 
detection efficiency of the leptons from the decay x*j ~^ ^ ^n- There is a factor 4 
appearing in the significances, since the asymmetries require two sets of equal beam 
polarizations V, as well as two decay modes, x°j ~~ * an d the charge conjugated 



decay x? 
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For an ideal detector, a significance of, e.g., S = 1 implies that the asymmetries 
can be measured at the statistical 68% confidence level. In order to predict the 
absolute values of confidence levels, clearly detailed Monte Carlo analysis including 
detector and background simulations with particle identification and reconstruction 
efficiencies would be required, which is however beyond the scope of the present 
work. 
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